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We show that if the hadronic resonance gas (HRG), with viscosity to entropy ratio n/s ~0.24, 
is physical at temperature T ~220 MeV, charged particles pr spectra and elliptic flow in Au+Au 
collisions at RHIC, over a wide range of collision centrality do not distinguish between initial QGP 
fluid and initial hadronic resonance gas. Unambiguous identification of bulk of the matter produced 
in Au+Au collisions require clear demonstration that HRG is unphysical at temperature T >200 
MeV. It calls for precise lattice simulations with realistic boundary conditions. 
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Experiments at Relativistic Heavy Ion Collider (RHIC) 
produced convincing evidences that in central and mid 
central Au+Au collisions, a hot, dense, strongly inter- 
acting matter is created p| , @ , @] , 0]. Whether the mat- 
ter can be identified as the lattice QCD predicted Quark 
Gluon Plasma(QGP) or not is still a question of debate. 
The problem is closely related to quark confinement; 
quarks are unobservable and any information about the 
initial state has to be obtained from observed hadrons 
only. QGP, even if produced in Au+Au collisions, is a 
transient state, it expands, cools, hadronises, cools fur- 
ther till interactions between the hadrons become too 
weak to continue the evolution. Hadronisation is a non- 
perturbative process. Whether or not the hadronisation 
process erases any memory of the constituent quarks is 
uncertain. If the hadronisation process erases the mem- 
ory, from the observed hadrons one can not comment on 
the initial QGP phase. Present search for QGP at RHIC 
is on the premise that the hadronisation process does not 
erase the memory and from the observed hadrons, using 
a dynamical model like hydrodynamics, one can back- 
trace to the initial QGP phase. Hydrodynamic equa- 
tions are closed with an Equation of State (EOS) and 
one can incorporate the possibility of phase transition 
in the model. Indeed, a host of experimental data in 
Au+Au collisions at RHIC is well explained in a hydro- 
dynamical model with QGP in the initial state Q. The 
alternative, namely hadronic resonance gas (HRG) in the 
initial state generally give poorer description to the data 
@, that too with HRG at very high temperature, Tj 
270 MeV. At such high temperature hadrons density is 
very large phad ~ 4 fm and it is difficult to believe that 
they retain their individual identity. However, situation 
changes if viscous effects are included. Model calcula- 
tions QiH^Bipil indicate that for a hadronic resonance 
gas, viscosity to entropy ratio is considerably larger than 
the ADS/CFT limit, r)/s > l/4n p). Since entropy is 
generated during evolution, unlike an 'ideal' HRG, a 'vis- 
cous' HRG can be initialized at a lower temperature such 
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that hadrons retain their identity yet reproduce the ex- 
perimental multiplicity. In the following, we show that if 
the hadrons retain their identity at temperature T «220 
MeV, charged particles pr spectra and elliptic flow in 
0-40% Au+Au collisions, do not distinguishes between 
initial HRG with viscosity to entropy ratio r//s «0.24 
and initial QGP with ?7/s=0.08. Unambiguous identifi- 
cation of the matter with QGP will require clear proof 
that HRG has a limiting temperature T < 200 MeV or 
the confincmcnt-dcconfinement transition occur at a tem- 
perature < 200 MeV. 

Most reliable information about the confinement- 
deconfinement transition temperature is obtained in lat- 
tice simulations of QCD. It is now established that the 
confinement-deconfinement transition is not a thermo- 
dynamic phase transition, rather a cross-over. Since it is 
cross-over, there is no unambiguous temperature where 
the transition take place. Inflection point of the Polyakov 
loop is generally quoted as the (pseudo) critical tempera- 
ture for the confinement-deconfinement transition. Cur- 
rently, there is debate over the value of the pseudo crit- 
ical temperature. HotQCD collaboration, with physical 
strange quark mass and somewhat larger than physical 
u and d quark masses {m s /m u _d = 10) claimed that 
both the chiral and confinement-deconfinement tran- 
sition takes place at a common temperature T c — 
192(4)(7). |l3,0,Q,[ii,[3. For more physical 
u and d masses (m s /m u . <j=20), the critical tempera- 
ture shift by ~ 5 MeV to lower side [T3|. Wuppertal- 
Budapest collaboration, with physical strange and light 
quark masses (m s /m u s = 28), on the other obtained 
a different result. [3, PI, (2Q1, (HI- Chiral transi- 
tion temperature T c =157 (3) (3) MeV is - 20 MeV less 
than the confinement-deconfinement transition tempera- 
ture T c =170(4)(4) MeV. Several other observables, e.g. 
quark susceptibility, pressure etc in Wuppertal-Budapest 
simulations is also shifted by 20-30 MeV compared to 
HotQCD. The ~ 20 MeV shift in the confinement- 
deconfinement transition temperature is not properly un- 
derstood. Possible reasons could be larger physical quark 
masses in HotQCD simulations. Also, one notes that 
though both the collaborations use staggered fermions, 
the fermion actions are different. It is still unproven 
that in the continuum limit, staggered fermions represent 
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QCD. The difference in the pseudo criticai temperature 
may also be due to different 'staggered' fcrmion action 
in the two simulations. Additionally, lattice results for 
the transition temperature can be up by ~ 30 MeV. All 
the lattice simulation uses periodic boundary condition, 
while in a realistic situation e.g. heavy ion collisions, 
the deconfined region is bordered by a confined phase. 
Exploratory quenched study suggests that with realistic 
boundary condition transition temperature can be up by 
30 MeV [22f . Considering the uncertainties associated 
with lattice simulations, the pseudo critical temperature 
for the connnement-deconfinement transition tempera- 
ture could be as high as T c ps 220 MeV. Qualitatively, 
for hadron size ps 0.5 fin, limiting hadron density (such 
that hadrons are not overlapped) is p\°^ it = 1/V ps 2 
/to -3 . For HRG, pi ia d ~ 2fm~ 3 corresponds to lim- 
iting temperature Tn m n= 220 MeV, value close to the 
highest possible pseudo critical temperature, as argued 
above. T «220 MeV is possibly the highest temperature 
at which hadrons can retain their individual identity. At 
higher temperature hadrons will overlap extensively and 
lose their identity. 

In the following we simulate Au+Au collisions in two 
scenarios, (i) no phase transition (NPT) scenario when a 
HRG with limiting temperature T=220 MeV is produced 
in the initial collisions and (ii) a phase transition (PT) 
scenario when QGP is produced. Space-time evolution 
of HRG/QGP fluid is obtained by solving, 

dpTi* = 0, (1) 
Dtt^ = - — -2r l W < > 1 u v> ) 

- [u"7r" A +u 1, ti vX ]Du x . (2) 

Eq[T] is the conservation equation for the energy- 
momentum tensor, T^ v = (e + p)u^u v - pg^ v + tt^ , 
£, p and u being the energy density, pressure and fluid 
velocity respectively, is the shear stress tensor (we 
have neglected bulk viscosity and heat conduction) . Eqf2] 
is the relaxation equation for the shear stress tensor tt^ v . 
In Eqf21 D = u^d^ is the convective time derivative, 
V<*V> = |(VV + V/) - \{d.u){g^ v - u^u v ) is 
a symmetric traceless tensor, rj is the shear viscosity 
and tv is the relaxation time. It may be mentioned 
that in a conformally symmetric fluid relaxation equa- 
tion can contain additional terms [23| . Assuming longi- 
tudinal boost-invariance, the equations are solved with 
the code 'AZHYDRO-KOLKATA' [13 , (H , (SI , (S3 in 
(t = y/t 2 — z 2 , x,y,r]= |^|^f ) coordinates. 

Hydrodynamic equations (EqQ] and [2]) are closed only 
with an equation of state (EOS) p — p(e). In the NPT 
scenario, we use EOS for the non-interacting hadronic 
resonance gas, comprising all the resonances with mass 
m res <2.5 GeV. In the PT scenario, we use a lattice 
based EOS [HI], where the high temperature phase is 
modeled by the recent lattice simulation [l3j]. At the 
cross-over temperature T co =196 MeV, the EOS smoothly 



changes to that of the hadronic resonance gas. 

Solution of Eq[T] and [5] require initial energy density, 
velocity distribution in the transverse plane at the initial 
time. A freeze-out prescription is also needed. In the 
PT scenario, we assume that QGP fluid is thermalised 
in the time scale Ti=0.6 fin 5]. Compared to QGP, a 
hadronic resonance gas is expected to thermalise late. 
We assume the canonical value, 7j=l fin as the initial 
time in the NPT scenario. Both in the PT and NPT 
scenario, initial fluid velocity is assumed to be zero. In an 
impact parameter b collision, the initial energy density 
is assumed to be distributed as 

e(b,a:,y) = eo[(l~f)N part (b,x,y) + fN coU (b 1 x,y)} (3) 

where N part (h, x, y) and N co u(h,x,y) are the transverse 
profile of participant numbers and binary collision num- 
bers. / in Eqf5]is the fraction of hard scattering. Most of 
the hydrodynamic simulations are performed with hard 
scattering fraction /=0.25 or 0.13 [5], 29]. We assumed 
/ = 0.25 in PT scenario. Hard scattering fraction is as- 
sumed to be zero in the NPT scenario. £o in Eqf5]is the 
central energy density of the fluid in impact parameter 
b=0 collision. In the PT scenario, for a fixed freeze- 
out temperature Tp=150 MeV, Eq is varied to best re- 
produce experimental charged particles pr spectra in 0- 
10% Au+Au collision. In the NPT scenario, we used the 
limiting value, £o=5.1 GeV/ fm 3 , corresponding to cen- 
tral temperature Ti=220 MeV. HRG initialised at cen- 
tral energy density £o=5-l GeV/ fm 3 do not produce ad- 
equate number of hadrons from the freeze-out surface 
at Tp=150 MeV. We lower the freeze-out temperature 
to Tp=110 MeV. Dissipative hydrodynamics also require 
initialisation of the shear stress tensors tt^ . In both 
the scenarios, we initialise ir^ at the boost-invariant 
values, n xx = tt vv = 2r](x,y)/3n, n xy = [27|. For 
the relaxation time T n , we use the Boltzmann estimate, 
t-k = 6r]/Ap. 

Results of our simulations for charged particles pt 
spectra and elliptic flow in 0-40% Au+Au collisions are 
shown in FigfT] and [D We have assumed that through 
out the evolution, viscosity to entropy ratio remains a 
constant. In PT scenario, we have simulated Au+Au col- 
lisions for four values of viscosity, (i)rj/s=0 (ideal fluid) 
(ii)j7/s = 1/4tt ps 0.08 (ADS/CFT limit), (iu)r)/s=0.12 
and (iv)?7/s=0.16. Corresponding central energy densi- 
ties are e =35.5, 29.1, 25.6 and 20.8 (GeV/fm 3 ) respec- 
tively. In the NPT scenario, we simulate Au+Au colli- 
sions for five values of viscosity, r//s=0, 0.08, 0.12, 0.24 
and 0.30. 

In Fig[TJ we have compared simulated charged par- 
ticles pt spectra in 0-40% Au+Au collisions with the 
PHENIX data [13]. The left panels (a-d) show the fit 
to the data in the NPT scenario. The solid lines from 
bottom to top corresponds to rj/s=0, 0.08, 0.16, 0.24 
and 0.30 respectively. Simulation results in the PT sce- 
nario are shown in the right panel (e-h). The lines from 
bottom to top are for viscosity to entropy ratio r]/s=0, 
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FIG. 1: PHENIX data [|3] for charged particles pr spectra 
in 0-40% Au+Au collisions are compared with hydrodynamic 
simulations in NPT (the left panel) and PT scenario (the right 
panel). See the text for details. 

0.08, 0.12 and 0.16 respectively. One observes that ideal 
HRG, initialised to central temperature T^=220 MeV, do 
not explain the pt spectra, data are largely under pre- 
dicted. Discrepancy between simulated spectra and ex- 
periment diminishes with increasing viscosity and data 
are best explained with rj/s=0.24:. Description to the 
data deteriorates if viscosity is further increased. In the 
PT scenario also, best fit to the 0-40% data is obtained in 
viscous QGP evolution, with r?/s=0.08-0.12. Ideal QGP 
or QGP fluid with rj/s > 0.12 give comparatively poorer 
description. To be quantative about the fit to the data 
in two different scenarios, we have computed x 2 values 
for the fits obtained to the 0-40% data, 

l E (e^)-^ 
TV ' erriiY 

i v 7 

We include both the statistical and systematic errors 
in the analysis. The X*/N values are shown in FigJSJi. 
In the NPT scenario, x 2 /TV as a function of r)/s shows 
a minima at 77/s=0.24, when (x 2 /-^0rmn=3.3. It is very 
interesting to note that the ?7/s=0.24, obtained from the 
analysis is in close agreement with theoretical estimates 
of viscosity to entropy ratio of a hot hadronic resonance 
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FIG. 2: (color online) same as in Fig[T]but for elliptic flow. 



gas, 7y/s=0.24-0.30 [7(. In the PT scenario also, best fit 
to the 0-40% data is obtained in viscous QGP evolution, 
with 7y/s=0.08, when (x 2 /N) min —5.6. Comparable fit 
is obtained with QGP viscosity ry/s=0.12. Qualitatively, 
with the initial conditions as used here, charged parti- 
cles px spectra in 0-40% Au+Au collisions are better 
explained with HRG than QGP in the initial state. 

In fig El we have compared simulated elliptic flow in 
0-40% Au+Au collisions with PHENIX measurements 
31]. PHENIX collaboration obtained V2 from two inde- 
pendent analysis, (i) event plane method from two inde- 
pendent subdetectors, v 2 {BBC} and v 2 {ZDC - SMD} 
and (ii) two particle cumulant w 2 {2}. All the three mea- 
surements of vi are shown. They agree within the sys- 
tematic errors. Simulated flows in the NPT scenario are 
shown in the left panels, the lines (top to bottom) cor- 
respond to HRG with 77/5=0,0.08, 0.16, 0.24 and 0.3 re- 
spectively. Flows in the PT scenario with QGP viscosity 
i]/s=0, 0.08, 0.12 and 0.16 are shown in the right pan- 
els. For elliptic flow also, we have computed \ 2 values 
for the fits. The results are shown in Fig(3jD. In the 
PT scenario, best fit to the 0-40% charged particles el- 
liptic flow data is obtained with QGP viscosity ?j/s=0.08 
((x 2 /TV) mi „=6.7). In the NPT scenario, the best fit is 
obtained with HRG with viscosity rj/s— 0.24, minimum 
(x 2 /TV) m j n =ll.l. Note that exactly, at these values of 
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FIG. 3: The filled circles and squares are x 2 /N for (a)pr 
spectra, (b) elliptic flow and (c) px spectra+elliptic flow, in 0- 
40% Au+Au collisions, in PT (initial QGP) and NPT (initial 
HRG) scenarios. 



viscosity, px spectra are also best explained. However, 
as opposed to the px spectra, elliptic flow is better ex- 
plained with QGP in the initial state rather than with 
initial HRG. If both the spectra and elliptic flows are 
analysed together, minimum y 2 is similar in both the 
scenarios, Xmin ~6-l in the PT scenario and «7.2 in 
the NPT scenario (see FigJ3t). It can not be claimed 
that experimental data are better explained with QGP in 
the initial state. Quantitatively, both the scenarios give 
nearly identical description to the data. Hydrodynamical 
evolution of initial QGP fluid with viscosity to entropy 
ratio ?7/s=0.08, thermalised in the time scale Ti— 0.6 fm 
to central energy density e =29.1 GeV/fm 3 and that of 
HRG with ?7/s=0.24, thermalised in the time scale Tj=l 
fm to central energy density £q=5.1 GeV/ fm 3 give nearly 
equivalent description to the PHENIX data for charged 
particles px spectra and elliptic flow in 0-40% Au+Au 
collisions. We have not explored all possible initial con- 
ditions in PT/NPT scenario. However, as the minimum 
X 2 /N value in NPT/PT scenarios is reasonably small, 
it is unlikely that with a different initial condition, data 
could be much better explained in the PT scenario, as 
long as HRG at central temperature 7^=220 MeV is phys- 
ical. 

One may argue that whether or not QGP is produced 
do not depend only on the px spectra and elliptic flow. 
Jet quenching, high px suppression etc. do give addi- 
tional support to the claim of QGP formation at RHIC 
energy collisions. However, high px suppression or jet 
quenching phenomena samples only a restricted phase 
space, as they are associated with high px trigger. Bulk 
of the matter are at low px- To claim that QGP is pro- 
duced in bulk, it is essential that charged particles px 
spectra and elliptic flow are explained in a hydrodynamic 
model. 

The present result that if HRG is physical at T=220 
MeV, NPT and PT scenario give equivalent description 




FIG. 4: (color online) In left panels (a)-(d) charged particles 
Pt spectra in 0-40% Au+Au collisions are compared with hy- 
drodynamical simulations with initial HRG and initial QGP. 
The solid and dashed lines corresponds to evolution of HRG 
with viscosity to entropy ratio r;/s=0.5, thermalised at r^=3 
fm, and QGP with viscosity to entropy ratio 77/3=0. 08, ther- 
malised as Ti=0.6 fm. In panels (e)-(h) results for elliptic flow 
are shown. 



to RHIC data on px spectra and elliptic flow will however 
be changed if hadronic resonance gas is physical only at 
a lower temperature T < 200 MeV. Even with large vis- 
cosity, r]/s=0.5, hydrodynamical evolution of HRG, ini- 
tialised at central temperature Ti=200 MeV, do not pro- 
duce comparable number of hadrons as in experiment, 
unless initial or the thermalisation time is large r» >3 
fm. In FigHl the solid lines are simulation results for 
charged particles px spectra and elliptic flow with vis- 
cous (t)/s=0.5) HRG with central temperature Ti=200 
MeV, thermalised at the time scale Tj= 3 fm. For com- 
parison, we have also shown in the simulations results 
with minimally viscous (77/5=0.08) QGP at the initial 
state (the dashed lines). One note that in 0-10% and 10- 
20% collisions, HRG with central temperature Ti=200 
MeV, thermalised in the time scale Ti=3 fm, give compa- 
rable description to the experimental charged particles 
px spectra, as obtained with QGP in the initial state. 
In 20-30% and 30-40% collisions however, data are bet- 
ter explained with QGP in the initial state. Similarly 
for the elliptic flow also. % 2 analysis also indicate that 
charged particles px spectra and elliptic flow in 0-40% 
Au+Au collisions are better expalined with initial QGP 
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than with initial HRG with central temperature Ti=200 
MeV. \ 2 /N «22 with HRG at the initial state, factor of 
~ 3 larger than that the value obtained with minimally 
viscous QGP in the initial state. Apparently, if HRG 
is physical only at temperature T <200 MeV, PHENIX 
data are better explained with initial QGP state than 
initial HRG state. Large thermalisation time Ti=3 fm, 
also raises the issue of non-equilibrium contribution to 
the particle production. Note that in dissipative hy- 
drodynamics, Tj is the time required by the system to 
achieve near equilibration. Large Tj imply that for a sub- 
stantial time, system remain in a highly non-equilibrium 
state, which can not be modeled theoretically, and contri- 
bution to particle production from this non-equilibrium 



state, though possibly can not be neglected, can not be 
estimated either. The highly non-equilibrium state, with 
unknown contribution to particle production make the 
HRG model with limiting temperature 7^=200 MeV, un- 
tenable. One can possibly exclude initial HRG if HRG 
is physical only at temperature T <200 MeV. Though, 
present lattice simulations indicate that HRG is physical 
only at <200 MeV, as discussed earlier, simulations are 
not precise enough to exclude HRG at T ss220 MeV. For 
confirmatory identification of matter produced in RHIC 
Au+Au collisions with QGP, it is essential to exclude the 
HRG scenario. Lattice simulations with realistic bound- 
ary conditions, preferably with Wilson fermion actions, 
are urgently needed to exclude possible HRG scenario. 
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